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ABSTRACT 


The objective of this study is to develop simplified micromechanical models to 
analyze the interface debonding between fiber and the matrix materials. Both analytical 
and simplified finite element models are used to predict the effective transverse elastic 
modules of fibrous composites with a partial interface crack based on the material 
properties of their constituents. The simplified finite element model uses springs in the 
connecting nodes between the fiber and matrix. A detailed finite element analysis, which 
is programmed using the MATLAB engineering software is performed to check the 
accuracy of the simplified models. The simplified models yield accurate effective 
transverse elastic modulii of various composites with partial interface cracks when 
compared to the results obtained from detailed finite element analyses. 
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1. INTRODUCTION 


Use of fibrous composite material is expanding rapidly, especially in structural 
application because properties of fibrous composite materials are beneficial such as high 
ratios of stiffhess-to-weight and strength-to-weight. 

Prediction of effective modulii of fibrous composites from material properties of 
their constituents has been a main topic, especially when there are micro-cracks. There 
are many different micromechanics models which have been proposed to predict effective 
moduli. However, most of them did deal with neither damages nor 3-D models. 
Broutman and Sahn observed the transverse matrix cracking in cross-ply composites 
[Ref. 1]. After the observation, the phenomenon has been studied anal 5 ftically and 
experimentally by several investigators. Highsmith and Reifsnider [Ref 2] conducted an 
experiment to measure the longitudinal stiffness reductions of various laminates. Talreja 
[Ref. 3,4] used continuum damage mechanics to represent damage as internal variables in 
a phenomenological approach, and Kwon [Ref 5] used an analytical micromechanics 
model to predict the effective modulii of composites. Kwon and Bemer [Ref. 6] also 
used a micromechanics model to investigate matrix cracking in laminated fibrous 
composite, and predict matrix cracking and its effect on stiffhess reduction. 

An analysis of a fibrous composite structure requires information of material 
properties of the composite. To determine the material properties of a composite, it is 
necessary to conduct experiments that will allow measurement of those properties. Even 
though some of the material properties can be determined with physical tests, it is 
difficult to obtain some other properties of the composite. 

This study will develop simplified micromechanical models to analyze the 
interface debonding between the fiber and the matrix materials using an analytical model 
and a simple two-dimensional finite element model with which the effective transverse 
elastic modulii of fibrous composites with partial interface cracks can be predicted. Since 
modeling an interface crack in the simplified finite element model is limited by a 
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minimum number of node points, the model uses springs in the connecting nodes 
between the fiber and matrix. In order to check the accuracy of the two simplified 
micromechanics models, a detailed finite element analysis is performed. 

Using properties of the fiber and matrix materials as well as their fiber volume 
fractions, it can be determined the relationship between the elastic modulus and the crack 
length. Figure 1 shows a micromechanical modeling process. In later chapters, both 
analytical and finite element micromechanical models will be presented, and followed by 
the results and discussion. The last chapter has conclusions and recommendations. 
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Figure 1. Modeling Process of a Micromechanical Model From a Fibrous 

Composite. 















II. MICROMECHANICS MODELS 


A. ANALYTICAL MICROMECHANICS MODEL 

This micromechanics model, which considers a unit cell consisting of the fiber 
and surrounding matrix material, was developed by Kwon [Ref. 5], and is shown in 
Figure 2. The cross-section of the fiber is assumed to be square, for simplicity, and the 
unit cell is divided into four subcells a,b,c and d. Subcell ‘n’ represents the fiber section 
and other subcells are the matrix. The size of the unit cell is unity, and the size of the 
subcell varies with the fiber volume fraction. The stress continuity at the subcell 
interfaces is expressed as: 

^22 ^22 ’ ^22 “ ^22 ’ ^33 ^33 ’ ^33 ” ^33 

^12 ” ^ 12’^12 "" ^12 ’ ^13 ”■ ^ 13’^13 ~ ^13 

a _ b _ a _ b 
^23 ““ ^23 ^23 ^23 ’ 

and the strain compatibility is assumed as: 


^11 “ ^11 ~ ^11 ^11 

S22 + ^22 ~ ^22 ^22 ’ ^33 ^33 


*33 ^*33 


£2 + £n = £n + £h > £n + 


The constitutive relation within each subcell can be expressed as: 


_cr __ T^a a 
^ if — ^ijkl^kl 5 


where i,j,k,l = 1,2,3 and a = a,b,c,d. 

From Figure 2 (before crack), Kwon [Ref 5]expressed the composite stresses and 
strains as a function of the subcell stresses, subcell strains and the fiber voltune fraction 


as follows: 
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where and f,^are composite stresses and strains, cr"y and f“;/ are stresses and strains 

of each subcell, and Vf is the fiber volume fraction that determines the size of the subcell. 
The composite stresses and strains are determined from the volume average of the 
stresses and strains of the subcells. 


2 

4 


d 

b 

c 

a 




-^r 


i-M 


a - Fiber Subcell 

b, c and d - Matrix Subcells 


Figure 2. Subcells of a Unit Cell. 


When there is an interface crack between the fiber and the matrix as shown in 
Figure 3, some volume portions of subcells a and c cannot sustain a load applied in the 3- 
direction. As a result, equation (2-4) is rewritten as 
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(2-6) 
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a - Fiber Subcell 

b, c and d - Matrix Subcells t Crack Length 


Figure 3. Crack Situation in Unit Cell. 


In order to find the relation between composite stresses and composite strains, 
first of all equations (2-1) through (2-3) and (2-5) are solved together. Then, the 
following matrix equation can be obtained. 

= = 1,2,3 , (2-7) 

where [A] is a constant coefficient matrix and is invertible. 

Pre-multiplying [A ]on both sides of equation 2-7, we get: 
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Finally, combining equation 2-4 and equation 2-8 along with equation 2-3 yields the 
relationship between composite stresses and strains: 

{5„}=[C]{s„} . (2-9) 

where [C] is the smeared material property matrix of the eomposite. In case of the 
interface crack, equation (2-6) is used instead of equation (2-4) with the same procedures. 

In Kwon's model, he assumed the fiber and matrix were orthotropic materials. 
This means the normal stress (strain) components and the shear stress (strain) components 
are uncoupled. Appendix A shows the equations relating the subcell strains to the 
composite strains for the normal strains. 

B. FINITE ELEMENT MICROMECHANICS MODEL 

The finite element method (FEM) is a numerical procedure for solving differential 
equations of physics and engineering. A continuous quantity such as displacements, 
temperature, or pressure can be approximated by a piece-wise continuous function 
spanned over a finite number of subdomains into which the problem domain is divided. 
For each subdomain, nodal variables are defined as unknowns which satisfy the 
governing equation in an approximate manner. The nodal unknowns are assembled into a 
system of matrix equation for the whole domain. The matrix equation is solved for the 
unknowns with prescribed boundary conditions. 

1. Finite Element Method 

In this study, one of the finite element methods called the weighted residual 
method (WRM) has been selected for the micromechanical model. The concept of this 
method is to minimize the residual multiplied by a weighting function (W). The WRM 
has many forms, and one that is very useful for solving boundary value problems is 
Galerkin's method. This is the method used in this study. The proeedure for this method 
is given below: 

The equilibrium equations for two dimensional elasticity are: 
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^<^x 


= 0 , 

dn 

dy 


dx 



yx 

+ 

= 0 , 

dx. 

4^ 



( 2 - 10 ) 


and the shear stresses satisfy: 


^ xy ^ yx » ( 2 “ 11 ) 

The two-dimensional strain-displacement relationships i.e. (kinemetic equations) are: 
du dv du dv 

= r,. =^ + — > (2-12) 


du 


_du 


where u and fi represent the displacements in x and y direction. 

The stress-strain relationships i.e. (constitutive equation) for the plane stress condition 
are: 


--^ V 1 0 

1 -// \-u 

0 0 y 

2 Vxy 


(2-13a) 


M=MW ■ 

Where E and // are elastic modulus and Poisson’s ratio. 

The traction boundary conditions are: 

<l>x = (^x^x + ^xy*iy 

<l>y = '^xynx^<^y^y ■ 

The basic formulation of the weighted residual method is: 

l=\RW^d^ = 0 (/ = U) , 

n 

where R is residual of the governing equation, and Wi is a weight function. 
By applying equation 2-10 to the weighted residual equation we get: 


(2-13b) 


(2-14) 


(2-15) 
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Using integrations by parts and combining the two I functions yields: 


- - + ^xv 

r 3c ^ ^ 

J 3(V, 3fV, 

3c ^ ^ 


m 1 




(2-16) 


(2-17) 


0 ^ 


n 0 ^ ^ \r^ 

dy dx. 


rn\m 


, (2-18) 


where (t)^ and are boxmdary tractions as stated in equation 2-14. Here, /2 denotes the 
domain and /’indicates the boundary. 

Plug gin g the strain-displacement equations into the constitutive equations and 
substituting the resultant equation into 2-18 yields: 


j 3c 
n 0 


3fV^ 31V, 
^ 3c 


1^1 “ It 

Z 1. 

dy 3c 




(2-19) 


In this model, linear triangular elements are used for the formulation. That is, each 
triangular element has three nodes and each node has two displacements u, v in the x and 
y direction so that each element has six degrees of freedom. The shape function that is 
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used for the triangular element is Hi (x,y). The displacements are interpolated using the 
shape functions. 


or in matrix form: 


/=I 


( 2 - 20 ) 


//, 0 0 //j 0 

0 //, 0 //, 0 //, 


( 2 - 21 ) 


The displacement vector is 2 x 1, cmd can be expressed as the shape function 
matrix (2x6) times the nodal displacement vector (6x1); as shown above. We define a 
matrix [B] that is obtained by taking first partial derivatives of the shape functions with 
respect to the coordinates, that is: 


dx. 


0 


0 

dH, 

0 


dx. 


dx 



0 

dH., 

0 



dy 




dH^ 



3c 

dy 

dx 


dx. 


( 2 - 22 ) 


Because the weighting function (Wi) of the Galerkin method is the same as the shape 
function: 




(2-23) 


the weighting fimction matrix in equation 2-19 can be written as: 
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or 

[A:,]{rf)=[iinD][fi]A(4 = {F4 , (2-27) 

where A is the area of the triangular element, t is the thickness of the element. For unit 
thickness, t is equal to 1. 

The right hand side of the equation 2-27 is called the element force vector {Fg} 
which is dependent on the surface traction boundary conditions and represents the total 
force at the nodal points in the x and y directions. Combining the element stiffiiess 
matrices into the system global matrix equation gives: 
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[KK] [dd] ={FF} 


( 2 - 28 ) 


The FEM used in this study is to calculate the elastic modulus of the composite 
material, whether there is a crack or not. From equation (2-28) we can solve for the 
nodals displacement {dd} by applying the displacement boundary condition. After we get 
tfie nodals displacements, stresses are computed for each element. The elastic modulus 
can be calculated from: 

. (2-29) 

where is the composite average stress and e ^ is the composite average strain. In this 
study, we applied a uniform displacement to the micromechanics model as shown in 
Figure 4 as the displacement boundary condition. 
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a - Fiber Subcell 

b, c and d - Matrix Subcells 


Figure 4. Unit Displacement Model. 

Two different kinds of finite element analyses are performed. The first one is a 
simplified model in which each subcell in Figure 4 is divided into two triangular 
elements. The other one is a detailed model for which each subcell is divided into many 
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niimber of finite elements. The latter model was used to check the accuracy of the 
simplified analytical and finite element models. 


2. Interface Crack 

When a partial interface crack between the fiber and the matrix exists, the crack 
caimot be modeled explicitly for the simplified finite element model because there is no 
node to represent the crack tip. As a result, it is assumed that there are elastic springs at 
the two comer nodes, which connect the fiber and the matrix. The double node technique 
is used for this purpose. Since the nodes have two degrees of freedom (x and y 
directions), the spring constants at every node are kx and ky. The stiffiiess matrix for the 
two dimensional spring element is; 


\ 0 -k, 0 

0 A, 0 ky 

-k. 0 A:, 0 


Q -k Q ky 


(2-30) 


By adding the stiffness matrix into the global system matrix, springs connect the 
fiber to the matrix. The spring constant determines the status of a partial interface crack. 
For example, if the spring constant is very large, there is no crack and if the spring 
constant is zero, there is a complete crack at the interface. An intermediate value of spring 
constant represents a partial crack. The proper spring constants need to be determined 
later for a given size of interface crack. 

C. MATLAB COMPUTER PROGRAM DEVELOPMENT 


The MATLAB engineering software is chosen to program the finite element 
procedure, and to calculate the composite material property of a composite with perfect 
bonding or partial cracks between the fiber and matrix. The following procedures are the 
outline of the MATLAB program development: 

1. Define the numbers of elements and nodes and input x, y coordinates of the nodes 
associated with each element using 'GRID'. 
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2. Input the material properties {E & v) of the fiber and matrix, and calculate the 
kinematic matrix [B], the material property matrix [D] and the element stiffness 
matrix [K^] using ‘TESTT’. 

3. Use ‘CHECK’ to assemble each element stiffness matrix [K^] into the system 
stiffness matrix [KK]. 

4. Use ‘MODTT’ to apply the geometric boundary conditions and to modify the 
system global stiffness matrix [KK]. 

5. Compute the displacement {dd} and stresses using “STRES.” 

6. From the main program “MAIN,” we obtain the composite material property, and 
compare the results by plotting. 

7. For the interface crack model, we add the spring stiffness matrix to the system 
stiffness matrix before step 5, and follow the remaining steps. 
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III. RESULTS AND DISCUSSION 


In this study, we considered three kinds of composites, Graphite/Epoxy, 
Glass/Epoxy and Boron/Epoxy. The properties of these composites are listed in Table 1. 
The relationship between the transverse elastic modulus {E) of each composite and the 
fiber volume fraction {Vf) for these three composites is shown in Figure 5. The fiber 
transverse elastic moduli of Glass and Boron are higher than that of Graphite, while the 
matrix transverse elastic moduli are the opposite. Therefore, when the volume fraction is 
less than 0.5, the transverse elastic modulus of Graphite/Epoxy is higher than the others, 
but when the volume fraction is greater than 0.5 the situation is reversed. Figure 6 shows 
the function that best fit the relation between transverse elastic modulus and the fiber 
volume flection in graphite/epoxy materials. The figures of the other materials are shown 
in Appendix C. (See Figure 36,Figure 37) 


Table 1. The Properties of Matrix and Fiber in the Graphite/Epoxy, Glass/Epoxy, 
and Boron/Epoxy Composites Materials._ 


Property 

Composite 

Graphite/Epoxy 

Giass/Epoxy 

Boron/Epoxy | 


Epoxy 

Matrix 

Glass Fiber 

Epoxy 

Matrix 

Boron Fiber 

Epoxy 

Matrix 

Transverse 

Elastic 

Modulus 

(Gpa) 

15 

5.35 

72.38 

2.75 

400 

2.75 

Poisson’s 

Ratio 

0.49 

0.354 

0.2 

0.35 

0.3 

0.35 


When the interface between the matrix and fiber begins cracking, the interface 
will be separated partially, as shown in Figure 7a. Figure 7b. In order to imderstand the 
stress distribution within the fiber when there is a partial interface crack, a uniform 
displacement was applied along the uncracked interface of the fiber as shown in Figure 
7b. Stress distributions for the three fibers with 50% or 60% crack length are shown in 
Figure 8 to Figure 13. The figures for other crack lengths are shown in Appendix D. (See 
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Figure 38 through Figure 43) From the figures, it is found that the stress distribution for a 
half of the crack length of the fiber down to the bottom is much lower when compared to 
other regions. The same results can be shown in the stress distribution over the both fiber 
and matrix with 50% of crack length ratio and 0.36 of fiber volume fraction for the three 
composites as seen in Figure 14 to Figure 16. The figures of other cases are shown in 
Appendix E, Figure 44 through Figure 49. Besides the part of the fiber we mentioned 
above, the stress distribution at the lower region of the left-hand side of the matrix is also 
negligible. Therefore, for the analytical model, the lower stress regions were removed as 
shown in Figure 17, and equation (2-6) was used to calculate the average composite 
stresses and transverse elastic modulus for a different crack length. 

The results of the analytical model and the detailed FEM model with Vf= 0.36 for 
the three composites are shown in Figure 18 to Figure 20. The results for other fiber 
volume fractions are shown in Appendix F, Figure 50 through Figure 55. They compare 
the reduced transverse elastic modulus ratio (Ec/Eq) for different crack length ratios and 
different fiber volume fractions. According to the results, when the crack length ratio 
increases, the transverse elastic modulus decreases, and the error between the FEM model 
and the analytical model for the Graphite/Epoxy composite is larger than those for the 
other two composites. This is a result of the stress ratio of the region which is neglected 
in the analytical model. The stress ratio at these regions for the Graphite/Epoxy 
composite is higher than those for the other two composites. So the error between the 
FEM and analytical models for the Graphite/Epoxy composite is somewhat larger. 
Nevertheless, these two models produce very compatible results and the analytical model 
can be used to predict the effective transverse elastic modulus for different fibrous 
composites with a partial interface crack length. 

As mentioned in Chapter II, we model springs in the connection nodes between 
the fiber and the matrix by adding the spring stiffness matrix into the system global 
matrix, as shown in Figure 21. For this model, the spring constants kx and ky are chosen 
from 5 X 10‘° to 10^. When the spring constant is equal to or larger than 5 x lO'®, it’s 
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assumed that there is no interface crack between the fiber and matrix, and when the 
spring constant is less than 10^, the complete debonding between the fiber and matrix 
would take place. Otherwise there is a partial debonding. 

Based on the above concept, we define two crack cases, low crack case and upper 
crack case, as shown in Figure 22. When there is no crack at point “A” and point “B” 
starts to crack, ka is equal to 5 x 10’° and kb changes from 5 x 10'° to 10^. This is the 
low crack case. On the other hand, when point “B” totally debonds, this means kb is equal 
to 10^, and ka changes from 5 x 10'° to 10^. This is the upper crack case. The results of 
these two cases for different materials with 0.36 of the fiber volume fraction are shown in 
Figure 23 through Figure 28. Appendix G, Figure 56 through Figure 67 shows the other 
fiber volume fractions. These figures illustrate that the spring constant and the crack 
length ratio form an exponential relationship. The best fitting function for each 
relationship is also shown in the figures. 

When the spring constant decreases, no matter what the case is, upper or lower, 
the partial crack ratio increases, and the elastic modulus will decrease. Figure 29 through 
Figure 31 show the relationship of the transverse elastic modulus ratio and the crack 
length ratio for the three models. The other results for different fiber volume fractions are 
shown in Appendix H, Figure 68 through Figure 73. According to these results, the 
simplified FEM micromechanical model with a proper spring constant can also predict 
the effective elastic modulus where there is a partial interface crack. 

In the above results, it is assumed the shape of the fiber is rectangular. If we 
change the shape of the fiber to be octagon instead, then we consider a new 
micromechanical model. Figure 32 shows the two different micromechanical models. 
We calculate the effective transverse elastic modulii of the new micromechanical model 
with the same material properties that used in the first simplified model, and compared 
the results with the first model. 

Table 2 shows the effective modulii (undamage) of these two simplified models 
for the three different composites in which fiber volume fraction is 0.36, and Figure 33 
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through Figure 35 shows the effective modulii of the new micromechanical model in 
cracking condition. 


Table 2. Effective Modulis (Undamage). 


Modulis 

Model 


Rectangular Fiber Shape 

Octagon Fiber Shape 

Transverse Elastic Modulus 
of Graphite/Epoxy 

7.638 X 10^ 

7.503 X 10^ 

Transverse Elastic Modulus 
of Glass/Epoxy 

6.098 X 10^ 

5.2147 X 10^ 

Transverse Elastic Modulus 
of Boron/Epoxy 

6.538 X 10^ 

5.5310 X 10^ 

Fiber Volume Fraction 

0.36 

0.365 ; 
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a - Fiber Subcell 

b, c and d - Matrix Subcells 



A Partial Crack Situation of Subcells 



(b) 


Fiber Part of Subcell in Cracking 


Figure 7. Partially Separated. 
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Stress 



Figure 8. Stresses Distribution Of Graphite Fiber With 50% Crack Ungth. 
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Stress 



Figure 9. Stresses Distribution Of Graphite Fiber With 60% Crack Length. 


25 



Stress 



Figure 10. Stresses Distribution Of Glass Fiber With 50% Craek Length. 
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Stress 



Figure 11. Stresses Distribution Of Glass Fiber With 60% Crack Length. 
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Stress 




Stress 



Figure 13. Stresses Distribution Of Boron Fiber With 60% Crack Length. 
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Stress 



Figure 14. Stresses Distribution Of Graphite/Epoxy With 50% Crack Length. ( Vf 

0.36) 
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Stress 



Figure 15. Stresses Distribution Of Glass/Epoxy With 50% Crack Length. (F/ - 

0.36) 
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Figure 16. Stresses Distribution Of Boron/Epoxy With 50% Crack Length. (Vf — 

0.36) 
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Transverse elastic modulus ratio: Ec/Eo 
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0.6 
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—": analy.model 
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Crack length ratio; C 


Figure 18. Transverse Elastic Modulis Ratio Of A Graphiteffipoiy Composite. (V> 

0.36) 
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Figure 23. Spring Constant Of A Graphite/Epoxy Composite With Vf = 0.36 (Low 

Crack Case). 
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Spring const ka 



Crack length ratio C 


Figure 24. Spring Constant Of A Graphite/Epoxy Composite With Vf = 0.36 (Upper 

Crack Case). 
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Spring constant kb 



Crack length ratio C 


Figure 25, Spring Constant Of A Glass/Epoxy Composite With Vf = 0.36 (Low 

Crack Case). 
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Spring const ka 



Figure 26. Spring Constant Of A Glass/Epoxy Composite With Vf = 0.36 (Upper 

Crack Case). 
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Spring constant kb 



Figure 27. Spring Constant Of A Boron/Epoxy Composite With Vf - 0.36 (Low 

Crack Case) 
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Spring const ka 



0.7 0.8 0.9 

Crack length ratio C 



Figure 28. Spring Constant Of A Boron/Epoxy Composite With Vf - 0.36 (Upper 

Crack Case). 



Transverse elastic modulus ratio Ec/Eo 



Figure 29.Transverse Elastic Modulis Of Graphic/Epoxy Composite {Vf -0.36). 
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Transverse elastic modulus ratio Ec/Eo 



Figure 31. Transverse Elastic Modulis Of A Boron/Epoxy Composite (V^ - 0.36). 
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Transverse elastic modulus ratio Ec/Eo 





Transverse elastic modulus ratio Ec/Eo 



Figure 35. Transverse Elastic Modulis Of Boron/Epoxy For New Model (Vj- = 

0.365). 
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IV. CONCLUSIONS AND RECOMMENDATIONS 


Both modified analytical and simplified FEM micromechanical models 
developed in this study can predict the effective elastic modulii of composite materials 
with a partial interface crack between the fiber and matrix. As a result of this study, the 
following conclusions and recommendations are made. 

A. CONCLUSIONS 


• A simple, analytical, micromechanical model that neglects the non-effective 
load-carrying zone, can predict the effective elastic modulii accurately. 

• A simplified finite element micromechanical model that contains interface 
springs can also predict the effective elastic modulii accurately. 

• Different fiber shape micromechanical models which have the same fiber 
volume fraction, and the same material properties, yield comparable results. 

• When the fiber elastic modulus is much larger than the matrix elastic modulus, 
the results of the simplified models are more accurate than those for the case 
where the fiber and matrix elastic modulii are not much different. 

B. RECOMMENDATIONS 

• The micromechanics models need to be coupled with an analysis computer 
program for general composite structures. 

• Non-elastic deformations of composites need to be included in the present 
models. 

• A follow-up research from this study would be to develop an algorithm for 
dete rmining spring constants, dependend on the state of the interface crack. 
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APPENDIX A. SUBCELL STRAINS EQUATIONS 
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♦Equations 10 through 12 <3 = yjv ^, for a subscript “a,b,c,d” are subcells 





APPENDIX B. MATRIX [B] 
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APPENDIX C. TRANSVERSE ELASTIC MODULUS AND BEST FITTING 

FUNCTION. 


X 10 




Transverse elastic modulus 



Figure 37. Transverse Elastic Modulus Boron/Epoxy and Best Fitting Function 
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Stress 


APPENDIX D. STRESSES DISTRIBUTION OF A FIBER. 



Figure 38. Stresses Distribution of a Graphite Fiber With 70% Crack Length 



Stress 



Figure 39. Stresses Distribution of a Graphite Fiber With 
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Stress 



Figure 40. Stresses Distribution of a Glass Fiber With 70% Crack Length 
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Stress 





Figure 42. Stresses Distribution of a Boron Fiber With 70% Crack Length 
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Stress 



Figure 43. Stresses Distribution of a Boron Fiber With 80% Crack Length 
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Stress 


APPENDIX E. STRESSES DISTRIBUTION OF A COMPOSITE. 



Figure 44. Stresses Distribution of a Graphite/Epoxy Composite With 50% Crack 

Length {Vj- 0.49) 
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Figure 45. Stresses Distribution of a Graphite/Epoxy Composite With 50% Crack 

Length (k/= 0.64) 
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Stress 



Figure 46. Stresses Distribution of a Glass/Epoxy Composite With 50% Crack 

Length (^}-= 0.49) 
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Stress 



Figure 47. Stresses Distribution of a Glass/Epoxy Composite With 50% Crack 

Length (f/= 0.64) 
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Stress 



Figure 48. Stresses Distribution of a Boron/Epoxy Composite With 50% Crack 

Length (f/= 0.49) 
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Stress 



Figure 49. Str esses Distribution of a Boron/Epoxy Composite With 50% Crack 

Length (Fy= 0.64) 
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APPENDIX F. TRANSVERSE ELASTIC MODULUS RATIO. 






Transverse elastic modulus ratio: Ec/Eo 



Figure 51. Transverse Elastic Modulus Ratio of a Graphite/Epoxy Composite With 

^^= 0.64 
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Transverse elastic modulus ratio: Ec/Eo 



Crack iength ratio : C 


Figure 52. Transverse Elastic Modulus Ratio of a Glass/Epoxy Composite With Vf 

=0.49 
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ulus ratio: Ec/Eo 



Crack length ratio : C 


Figure 53. Transverse Elastic Modulus Ratio of a Glass/Epoxy Composite With 

=0.64 
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modulus ratio: Ec/Eo 







Transverse elastic modulus ratio; Ec/Eo 



Crack length ratio: C 


Figure 55. Transverse Elastic Modulus Ratio of a Boron/Epoxy Composite With Vf 

= 0.64 
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Spring constant kb 


APPENDIX G. SPRING CONSTANT OF A COMPOSITE. 



Figure 56. Spring Constant of a Graphite /Epoxy Composite with I^=0.49 (Low 

crack case) 
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Spring constant kb 



Figure 58. Spring Constant of a Glass /Epoxy Composite with Fjr-0.49 (Low crack 

case) 
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Spring const ka 



Crack length ratio C 


Figure 62. Spring Constant of a Graphite /Epoxy Composite with i^=0.49 (Upper 

crack case) 
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Spring const ka 



Figure 63. Spring Constant of a Graphite /Epoxy Composite with i^=0.64 (Upper 

crack case) 
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Spring const ka 



Crack length ratio C 


Figure 64. Spring Constant of a Glass /Epoxy Composite with 1^=0.49 (Upper crack 

case) 


87 







Spring const ka 



Figure 66. Spring Constant of a Boron /Epoxy Composite with f^=0.49 (Upper 

crack case) 
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Spring const ka 
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\ 


o \ 


\ 

O' 






Q _ 


O. _ 


o _ 




~ -o 


0.6 


0.7 0.8 0.9 

Crack length ratio C 


1.1 


Figure 67. Spring Constant of a Boron /Epoxy Composite with J^-0.64 (Upper 

crack case) 
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Transverse elastic modulus ratio Ec/Eo 


APPENDIX H. TRANSVERSE ELASTIC MODULUS OF COMPOSITE. 



Crack length ratio C 


Figure 68. Transverse Elastic Modulus of a Graphite/Epoxy Composite of Three 

Models {Vf = 0.49) 
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Figure 70. Transverse Elastic Modulus of a Glass/Epoxy Composite of Three Models 

{Vf = 0.49) 






Transverse elastic modulus ratio Ec/Eo 



Figure 71. Transverse Elastic Modulus of a Glass/Epoxy Composite of Three Models 

(Vf = 0.64) 


94 





Transverse elastic modulus ratio Ec/Eo 



Figure 72. Transverse Elastic Modulus of a Boron/Epoxy Composite of Three 

Models {Vf = 0.49) 
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Figure 73. Transverse Elastic Modulus of a Boron/Epoxy Composite of Three 

Models (F. =0.64) 
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